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Abstract 

We investigate the accuracy of first-principles many-body theories at the nanoscale by comparing 
the low energy excitations of the carbon fullerenes C20J C24, C50, Ceo, C70, and Cso with exper- 
iment. Properties are calculated via the GW-Bethe-Salpeter Equation (GW-BSE) and diffusion 
Quantum Monte Carlo (QMC) methods. We critically compare these theories and assess their 
accuracy against available photoabsorption and photoelectron spectroscopy data. The first ioniza- 
tion potentials are consistently well reproduced and are similar for all the fullerenes and methods 
studied. The electron affinities and first triplet excitation energies show substantial method and ge- 
ometry dependence. These results establish the validity of many-body theories as viable alternative 
to density-functional theory in describing electronic properties of confined carbon nanostructures. 
We find a correlation between energy gap and stability of fullerenes. We also find that the elec- 
tron affinity of fullerenes is very high and size- independent, which explains their tendency to form 
compounds with electron-donor cations. 

PACS numbers: 78.66.Tr,73.22.-f,78.20.Bh,78.40.Ri 
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INTRODUCTION 



Carbon fullerenes form a remarkable series of molecules whose physical and chemical 
properties are relatively well known. Their properties can be controlled through size, shape, 
and structure if]. They can be functionalized, doped with endohedral or exohedral atoms, or 
polymerized 3j. Ceo is the most stable and well-studied molecule in the series. It has been 
widely used e.g. in the design of efficient organic photovoltaic cells and as precursor in 
the growth of carbon nanotubes and nanopeapods 6|, [7|. The strength of carbon sp 2 bonds 
makes fullerenes extremely stable, and the ability of carbon to form bonds with different 
hybridization levels makes possible the existence of both hexagonal and pentagonal rings 
and hence a rich variety of sizes and shapes [l, 3]. In fact, carbon is the main constituent of 
a wide range of nanostructures from diamondoids to organic molecules, including fullerenes 
and nanotubes. It is likely that nano-electronic devices with desired properties can be 
created through intentional selection of specific carbon nanostructures [8|. In order to 
systematically design such devices it is necessary to develop a thorough understanding of the 
electronic and optical properties of the nanostructure, either isolated or in a functionalized 
medium. Detailed knowledge of the electronic structure of fullerene molecules is the key to 
explore several properties of those systems. As an example, the conductance of fullerene- 
based junctions depends crucially on the relative energy of electronic orbitals on the molecule 
with respect to the chemical potential on the leads |9| . 

Experimental characterization is often hindered at the nanoscale and therefore theory 
plays an increasingly important role in nano-science as a guide for experimental interpreta- 
tion. While it is often claimed that theory can be used to design nanomaterials with targeted 
properties, for theory to have predictive power it is crucial to validate it in those precious 
cases where the comparison with experiment can be made. Carbon compounds including 



nanotubes, 



10 



llj adamatanes and fullerenes are among the few nanomaterials that have 



been characterized in terms of structure and composition. 

In this paper we apply the two most accurate approaches for calcula ting electronic/optical 
properties, the GW-Bethe Salpeter equation (GW-BSE) approach 12, Q and diffusion 
quantum Monte Carlo (QMC) 3]. These methods are complementary: the GW-BSE ap- 
proach is based primarily on approximate Green's functions and in-principle yields the entire 
optical spectrum, while QMC uses a stochastic sampling approach that is in-principle exact 
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for the lowest energy excitation of each symmetry 14j . Due to recent algorithmic advances 
and a favourable scaling with system size, both of these methods can provide the optical 
properties of nanostructures with one hundred atoms or more using existing parallel com- 



puters 
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molecules 



12 



16]. While GW-BSE and QMC have previously shown good agreement for small 
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171 ]. this is the first comprehensive study of these methods for these 



much larger yet prototypical nanostructures. Owing to the complexity involved in both 
QMC and GW-BSE when applied to electronically confined systems, not much is known 
about the accuracy of either theory in confined nanostructures. 

We study a range of fullerenes ( C20, C24, C50, Ceo, C70, and Cso ), with the aim of 
answering these questions: (1) Do these methods accurately predict the optical and electronic 
properties, particularly for the best experimentally characterized Cqq molecule? (2) What are 
the trends in the physical properties across the fullerene series? (3) Where are improvements 
in each method required? 

METHOD 



Fullerenes are routinely produced in a wide range of sizes, from C20 to C100 or larger, 
and many isomers have been identified 18|. The existence of several low-energy isomers 
for each molecular size poses a difficulty in modeling fullerenes. For each computed size we 
have chosen high-symmetry structures that satisfy the isolated pentagon rule [l| as closely 
as possible. Figure [1] depicts the isomers investigated here, with the corresponding point 
groups. The atomic coordinates were obtained by minimizing the total energy within norm- 
conserving pseudopotential calculations using DFT and the Perdew-Burke-Ernzerhof (PBE) 



functional 



19(. Choices of grid point separation, energy cutoffs and supercell sizes were 



such that Kohn-Sham eigenvalues of bound orbitals are converged to 0.02 eV or better and 
bond lengths to 1% or better. Atom coordinates are included in the supporting information 
accompanying this paper ^ . 

Owing to its simplicity, density-functional theory (DFT) has been used to predict the 
electronic properties of nanostructures. Nevertheless, it has two major deficiencies: (1) It can 
only describe the energy of optical excitations through ad hoc approaches such as constrained 



DFT 



201 ]: and (2) it is not possible to systematically improve the description of exchange 



and correlation effects within DFT because the exact exchange-correlation functional cannot 
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be calculated in closed form [20j. In contrast, many-body approaches can be systematically 
improved, at the expense of increased numerical complexity 



In the framework of the GW approximation 13], electron self-energies are calculated in 



terms of the Green's function G and the screened Coulomb interaction W as: 

£(l,2) = iG(l,2)W (l,2) , (1) 

where screening is commonly evaluated in the random-phase approximation (RPA). We 
can modify Eq. (00) by replacing the RPA with time-dependent DFT with an adiabatic, 
local-density kernel (TDLDA). In order to retain Feynman diagrams with the same order, 
we add a properly symmetrized vertex term together with the TDLDA screened Coulomb 



interaction 



1$: 



S(l,2) =zy"d(3)G(l,3)W / (3,2)r(l,2;3) + symm. . (2) 

At its lowest level of approximation, the Green's function and the screened Coulomb inter- 
action above are coated within DFT Q. Further improvement in the theory is attained 
by imposing self-consistency among the self-energy E, Green's function G and screened 
Coulomb interaction W [3] We impose self-consistency in an iterative scheme. In the first 
iteration, we calculate the self-energy using DFT electronic orbitals and Eq. (T5]). In subse- 
quent iterations, we use a scissors operator to correct the energy of DFT orbitals, according 
to the self-energy calculated in the previous iteration. We stress the fact that the scissors 
operator is constructed using information from the previous iteration only. No free param- 
eters are involved. We stop when electron affinity and ionization potential are converged 
to within 0.02 eV. Usually, this is reached in less than five iterations. Tests performed on 
smaller molecules ( CI2, benzene, CO, SiH-4) indicate that going behond the scissors opera- 
tor method would improve the final results by 0.1 eV or less. Scissors operators have been 
used extensively in the context of electronic band structure of bulk semiconductors (see Ref. 



13 



and references therein). In the following, we refer to Eqs. ([I]) and ([2]) as GWo and 
GWj approximations respectively. We denote by scGW/ approximation the self-consistent 
solution of Eq. ([2]), where self-consistency is imposed through a scissors operator. After 
the quasi-particle orbitals are determined, we calculate optical excitations by diagonaliz- 
ing the Bethe-Salpeter equation for electrons and holes following standard methodologies 
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12L Il5l . Il7j |. The BSE method has been used successfully to predict excitation energies of 
organic molecules 15], CdSe nanocrystals 21] and hydrogenated silicon clusters 

In the QMC methodology (for a review, see Ref. 14J), excitation energies are obtained 
via the difference in total energy between the ground state and individually constructed ex- 



cited states. We use the diffusion Monte Carlo (DMC) technique with trial wave-functions 



consisting of a single S 
correlation function 



ater determinant of single-particle orbitals multiplied by a Jastrow 
221 ] . The most significant approximation is the fixed-node approxima- 
tion, whereby the Fermion nodes of a trial wave-function are used in place of the exact 
nodes, introducing a variational error. Although DMC calculations for large systems addi- 
tionally require the use of pseudopotentials, for carbon-based systems this approximation 
is secondary. We evaluate the pseudopotentials within the localit y ap proximation. A test 
calculation for C20 using a variational pseudopotential evaluation |23j found a global shift 
of 0.8 ± 0.1 eV but identical energy differences to within ±0.1 eV, demonstrating that this 
error is small in the fullerenes. In practice, the accuracy of the results obtained for excita- 
tions is solely determined by the accuracy of the ground and excited state nodal surfaces. 
In contrast, the GW calculations are dependent on the detailed shape of the Kohn-Sham 
orbitals and the initial Kohn-Sham eigenvalues. Given exact nodes for the ground state and 
for the first excited state of a given symmetry, DMC yields in-principle exact results In 



241 ] . In our calculations we use 



general, there are additional symmetry related qualifications 
the Kohn-Sham orbitals from the ground state DFT calculation in the Slater determinant 
to effectively define the nodal surface of the trial wavefunction. DMC was used recently 
to calculate excitation energies and ionization energies of diamondoids yielding very good 
results {25)]. 

To make comparisons between GW-BSE and QMC on an equal footing we used the same 
norm conserving pseudopotential in each set of calculations. The initial QMC wave-functions 
were obtained from a local-density approximation (LDA) DFT calculation, while for GW 
calculations we additionally investigated the PBE 19(] gradient corrected functionals. GW/- 
electron affinities and vertical ionization potentials obtained using the LDA are higher than 
the ones obtained using PBE by a fraction of electron-volt: 0.1 eV in Cqq and 0.3 eV in C2o- 
The impact of replacing the LDA with PBE on Kohn-Sham eigenvalues is even smaller: less 
than 0.05 eV. In view of these small differences we concentrate on LDA-derived results. 
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RESULTS 



We have investigated two isomers of Cso- DFT calculations have already indicated that 
Cgo exists in a large number of isomers, where the oval shaped structure (point group D 5( j) 



is the most stable one 



26 



27| . Our DFT calculations predict that the icosahedral isomer 



(lh) is more energetic than the Dsd isomer by 0.8 eV, confirming earlier calculations 26]. 
This difference increases to 1.74 ± 0.10 eV in QMC. 

Table [J shows the vertical ionization potential (IP) calculated within the GW and QMC 
theories. The ionization potential was also calculated using DFT only, as the difference in 
total energy between the neutral and the charged molecules ("ASCF-DFT") column. The 
IP is approximately 7.5 eV for all of the theories, with no well defined trend in IP with size 
apart from a slow decrease with increasing number of atoms. The analyzed structures have 
different curvatures, caused by the different positions of pentagonal rings from molecule to 
molecule. Therefore, geometry and symmetry become important factors in determining the 
ionization potential. A classical electrostatic model based simply on a c harg ed hollow sphere 



60 



can only predict qualitatively the ionization potential of fullerenes [28|, [29( . Except for G 
and C70, the QMC ionization potential is lower than the GW/ ionization potential. The 



discrepancy between the two approaches is around 0.2 eV, close to the numerical accuracy 



associated with each theory 



44[ | . Except for GWq all the methods predict IPs close to 



experiment. Self-consistency in the GW method is found to give very small improvement 
on the already accurate IP. Table [I] also shows that the GW/ approximation predicts an IP 
much more accurate than the GWo approximation. Therefore, including vertex corrections 
and polarizability effects beyond RPA are an improvement on the theory, confirming similar 
observations in benzene and naphthalene molecules 151 ] . 



Ultra-violet photoelectron spectroscopy (UPS) can determine not only the first ioniza- 
tion potential but also higher ionization energies, thus givin g di rect access to the electronic 
structure of the material. Lichtenberger and collaborators |30(] have identified five ioniza- 
tion bands in the Cgo spectrum. Table [IT] shows the assignment of the ionization bands 
and comparison between measured and calculated vertical ionization energies. Higher-order 
ionization energies can only be calculated within DFT if a minimization constraint is used. 
Calculation results presented on column "CDFT" on Table [III were obtained by minimizing 
the total energy of the cation Cg with the constraint that the ejected electron occupied a 
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TABLE I: Vertical ionization potentials calculated from DFT, QMC and GW theories. Error bars 
in QMC values are indicated between parentheses. All energies in eV. 





ASCF-DFT 


QMC 


GW 


GW f 


scGW/ 


Exp. 


C20 


7.31 


7.27(11) 


7.99 


7.35 


7.41 




C24 


7.77 


7.70(10) 


8.49 


7.86 


7.81 




C50 


7.29 


7.29(14) 


7.97 


7.33 


7.35 


7.61 a 


C60 


7.61 


7.86(21) 


8.22 


7.70 


7.86 


7.6 6 


C70 


7.54 


7.69(12) 


8.12 


7.53 


7.45 


7.47 c 


Cso (D5d) 


6.67 


6.30(10) 


7.24 


6.59 


6.65 


6.84 a 


Cso 0-h) 


6.86 


6.91(10) 


7.45 


6.90 


6.95 




average error 


-0.10 


-0.09 


0.51 


-0.09 


-0.05 




root mean square error 


0.18 


0.36 


0.52 


0.20 


0.21 



















a Ref. 


18 


. b Ref. 


30 


c Ref. 


rr> 1 





specific molecular orbital (HOMO, HOMO-1, or deeper orbitals). Whereas constrained DFT 
predicts very accurately the first two ionization energies, it increasingly underestimates the 
higher ionization energies. The discrepancy seems to be more than 1 eV in the ionization 
band at 12.4 to 13.8 eV. In contrast, the self-consistent GW/ approximation predicts ioniza- 
tion energies no more than 0.35 eV away from the experimental data over the same energy 
range. 

The electron affinity (EA), shown in Table HTT1 is also predicted differently by QMC and 
the various levels of GW theory. The former theory underestimates the electron affinity by 
as much as 1 eV relative to the latter one. Inclusion of vertex corrections (from GWo to 
GWf) reduces the electron affinity by approximately 0.6 eV. Inclusion of self-consistency 
(from GWf to scGW/) increases the self-energy by a further 0.5 eV, bringing the calculated 
results in very good agreement with experimental data and with ASCF-DFT predictions. 
The electron affinity is found to be around 2 to 4 eV. All theories employed in this work 
predict that C20 has the smallest electron affinity whereas CsoO-h) has the largest electron 
affinity. C^q is well known to be easily ionizable 1, jj, 31, 32|. The present results also show 
that less stable fullerenes are easily ionizable as well, which opens the path to functionalize 
fullerenes of different sizes. 
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TABLE II: Orbital assignment and vertical ionization energy for the highest occupied molecular 
orbitals in Cgo- Column "C-DFT" shows constrained DFT results. Columns GWo and GW/ are 
non-self-consistent GW calculations, with and without vertex corrections respectivelly. Column 
scGW/- includes self-consistency and vertex corrections. Experimental data quoted from Ref. I3C . 



Experimental data spanning multiple rows indicate ionization bands that could not be resolved. 
All energies in eV. 





C-DFT 


GW 


GW/ 


scGW/ 


Exp. 


h u 


7.61 


8.22 


7.70 


7.86 


7.6 




8.78 


9.33 


8.85 


9.23 














8.95 




8.90 


9.42 


8.93 


9.25 




K 


10.47 


11.93 


11.40 


11.93 






10.50 


11.00 


10.51 


11.07 


10.82-11.59 




11.03 


11.46 


11.02 


11.6 






10.79 


12.19 


11.65 


12.2 






11.79 


13.23 


12.73 


13.36 














12.43-13.82 


ti 9 


12.28 


13.60 


13.11 


13.76 






12.33 


12.74 


12.32 


13.04 





One important result of this work is that GWo an d GW j are shown to be consistent in the 
"electronic gap" , denned as the difference between ionization potential and electron affinity. 
This definition of electronic gap is the closest analogue of band gap to isolated molecules 
\\\ . Although GWo an d GW/ predict different values for IP and EA, the differences cancel 
out in the electronic gap. The reason for this cancellation is that contributions arising 
from LDA screening, such as the vertex T, have similar magnitude when evaluated at the 
highest occupied molecular orbital (HOMO) and at the lowest unoccupied molecular orbital 
(LUMO). As a consequence of this cancellation, the contributions from LDA screening are 



often neglected in calculations of electronic gap in crystals 13]. In molecular systems, these 



contributions are not negligible because they affect IP and EA separately 15|, |33[. Self- 



consistency has been found to systematically increase the electronic gap compared to non 



self-consistent GW methodologies in a variety of different materials 13, |34j . Our calculations 



confirm the same behavior in finite systems, and they also indicate that self-consistency 
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TABLE III: Electron affinity, or vertical detachment energy, calculated from DFT, QMC and GW 
theories. Error bars in QMC values are indicated between parentheses. All energies in eV. 





ASCF-DFT 


QMC 


GW 


GW/ 


scGW/ 


Exp. 


C20 


2.17 


1.76(11) 


3.55 


2.92 


2.36 


2.25 a 


C24 


3.04 


2.57(11) 


4.19 


3.55 


2.88 




C50 


3.73 


3.52(14) 


4.75 


4.12 


3.73 


< 3.10 6 


C60 


2.94 


2.23(19) 


3.87 


3.33 


2.98 


2.69° 


C70 


2.96 


2.46(11) 


3.98 


3.35 


2.83 


2.76 d 


Cso (h>5d) 


3.46 


3.25(10) 


4.62 


3.91 


3.88 


3.70 6 


Cso 0-h) 


3.98 


3.90(11) 


5.17 


4.61 


4.38 




average error 


0.03 


-0.42 


1.25 


0.53 


0.16 




root mean square error 


0.20 


0.43 


1.16 


0.56 


0.18 























a Ref. 


39 


. 6 Ref. 


40 


. c Ref. 


31 


. d Ref. 


32 





affects the energy of unoccupied orbitals more than of occupied orbitals. Self-consistency 



15 



21 



13, 



could be more important in fullerenes than in most small molecules studied so far 
I35I because of the narrower HOMO-LUMO gap. 

There is an ongoing debate in the literature concerning the importance of vertex correc 



tions 



The GW/ approximation has vertex corrections included in two ways: 
explicitly through the function T in Eq. ([2]), and implicitly through the screened Coulomb 
interaction (see e.g. Eq. 17 and 18 in Ref. 15|). We have observed that most of the dif- 
ference between GW/ and GWo calculation results originates from the explicit contribution 
rather than the implicit contribution. Matrix elements of the "explicit vertex contribution" , 
defined as AE = i(GWT — GW) have approximately the same magnitude at the HOMO 
and the LUMO of each molecule. As a function of molecule size, AE fluctuates around 
0.6 eV to 0.7 eV across the family of fullerenes studied, with the exception of C20 where it 
is around 0.8 eV. The "implicit vertex contribution", i(GW — GWo) is typically five times 
smaller in magnitude. Self-consistency following the prescription presented here does not 



modify significantly the strength of vertex contributions. T 
observed in small oligoacenes (benzene and naphthalene) 



lis behavior is similar to the one 



151 ] but it is somewhat different 



from the scenario observed in silicon nanocrystals. There, the explicit vertex contribution 
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TABLE IV: Excitation energy of the first spin-triplet state. Calculated values under the "GWj- 
BSE" column were obtained using the GWf approximation in the electron self-energy. Column 
"scGW/-BSE" indicates results obtained with the self-consistent GW/ approximation. Error bar 
in QMC values are indicated between parentheses. All energies in eV. 



TDLDA 



QMC 



GW/-BSE 



scGW/-BSE 



Exp. 



C20 


0.57 


0.87 (8) 


-0.12 


0.62 




C24 


0.43 


0.90(14) 


-0.07 


0.67 




C50 


0.20 


0.61(13) 


-0.17 


0.24 




C60 


1.52 


2.34(19) 


1.14 


1.65 


1.58° 


C70 


1.62 


2.31(12) 


1.01 


1.45 


1.55 6 


Cso (D5d) 


0.41 


0.07(10) 


-0.15 


0.20 




Cso 0-h) 


0.11 


0.26(10) 


-0.35 


-0.06 





average error 



0.01 



0.76 



-0.49 



-0.02 



root mean square error 


0.05 


0/ 


m 




0.49 


0.09 


a Ref. 


11 


. b Ref. 


12 





is no more than around 0.3 eV and cancellation between explicit and implicit contributions 
is more complete. 

Table IIVI shows the calculated and measured excitation energy of the first spin-triplet 
electronic states in each fullerene. Experimental data were obtained from phosphorescence 
decay measurements, and they include a Stokes shift absent in the theoretical calculations. 
Based on ASCF-DFT calculations, we estimate the Stokes shift to be no more than 0.2 eV. 
Independent evidence for small Stokes shift is found in the photoelectron spectra of C^~ 
and C^, which show a high 0-0 line followed by lower vibrational side bands at higher 



energy 



31 



32|. TDLDA gives good excitation energies, with contrasts with the well-known 



gap underestimation of carbon nanotubes 



361 ] and of periodic systems in general [12|. Self- 



consistent GW-BSE is also very accurate, having a discrepancy from experimental data 
comparable to the Stokes shift. The QMC data is too high by approximately 0.8 eV, and 
the non-self-consistent GW-BSE data is too low by approximately 0.5 eV. 

The difference in the QMC results must result from a poor cancellation of nodal error 
between ground and triplet excited states. For small molecules it is possible to improve 
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the nodal surface by constructin g m ultideterminant trial wavefunctions or performing addi- 



tional wavefunction optimization 37|, but these techniques are not yet practical for the large 
fullerenes. Clearly fullerenes represent a good benchmark for methods to improve the nodal 
surfaces. 

The GW-BSE results reflect the discrepancies in the electron affinity: since low-energy 
optical excitations usually involve the LUMO, excitation energies obtained without self- 
consistency are underestimated by approximately 0.6 eV (column "GW/-BSE" in Table 
E}. We further examine the observed discrepancies by computing higher-energy neutral 
excitations of Ceo- Table [V] presents excitation energies of the first few spin-singlet states 
and the projection of each excited state on the subspace of electron transitions from the 
HOMO to the LUMO. This projection is simply a sum over overlap integrals between the 
GW-BSE eigenstates |iS) and the many-body states obtained by promoting one electron 
from the HOMO to the LUMO in the ground state \G): 

HOMO LUMO 

p H- L = E E \{s\4^\g)\ 2 , (3) 

j i 

where a] and cij are many-body creation and destruction operators respectively. This pro- 
jection is not readily available in DMC since only the probability density consisting of the 
product of the QMC ground state and DFT-based trial wavefunction is available. 

Identification of the various excitations in Cqo is facilitated because the molecule is 
highly symmetric, leading to high degeneracy and wide energy separation between states. 
The HOMO and the LUMO of Cqo have symmetry representation h u and t* respectively. 
Electronic transitions between them give rise to four multiplets T*, H 3 , T^, G g . Table |V] 
shows that the excitation energy of all four multiplets is underestimated within the non- 
self-consistent GW/ approximation, typically by 0.6 eV. This is a direct consequence of the 
severe overestimation of the electron affinity if self-consistency is not imposed. Some of the 
higher-energy excitations do not involve the LUMO, such as excitation T^. For them, the 
GW-BSE prediction is compatible with experimental data. 

Finally, Table HVl shows that the most stable fullerenes, namely C 60 (1^) and C 70 (D5/J, 
have the largest triplet excitation energies. Less stable fullerenes have small triplet excitation 
energy, of the order of 0.4 eV, whereas CgoO-h) is predicted to have the smallest triplet energy: 
0.11 eV (TDLDA). Knowing that triplet excitation energies are related to the energy gap, 
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TABLE V: First spin-singlet excitation energies in Cgo- All energies in eV. We use the same 
notation of Table IIVI Column labeled "H-L" indicates the projection of this excitation onto the 
HOMO-LUMO subspace (see text). Experimental data from Ref. |43. 





TDLDA 


GW/-BSE 


scGW/-BSE 


H-L 


Exp. 


T 1 

9 


1.69 


1.47 


1.88 


99 % 


1.82-1.97 


rp2 

9 


1.71 


1.53 


1.92 


99 % 


1.82-1.97 




1.64 


1.41 


1.82 


99 % 


2.00-2.21 


H 9 


1.79 


1.70 


2.15 


94 % 


2.25-2.32 



this trend shows perfectly correlation between stability and energy gap: stable fullerenes 
have wide gap between HOMO and LUMO, unstable fullerenes have narrow gap between 
HOMO and LUMO. 



CONCLUSION 



We have calculated and analyzed the low energy excitations of fullerenes C20 to Cgo using 
three theoretical methods: density functional theory, many-body GW-BSE and diffusion 
Quantum Monte Carlo. Overall these methods give an accurate description of the opti- 
cal properties. Ionization potentials are approximately constant and well predicted by all 
methods. The more stable fullerenes (Cgo and C70) are found to have highest excitation 
energy. Stability correlates well with energy gap. We find a weak dependence of electron 
affinity with molecule size, with important implications for doping and functionalization of 
fullerenes. We find that two ingredients are essential to bring GW-BSE predictions to agree 
with experimental data: vertex corrections, which are responsible for an almost rigid shift 
of electronic orbitals with respect to the vacuum level, and self-consistency, which widens 
the energy gap between occupied orbitals and unoccupied orbitals. QMC calculations of 
the lowest energy triplet excitation using a single determinant of Kohn-Sham orbitals give 
results higher than experiment by approximately 0.8 eV, indicating that the DFT-derived 
trial wave-functions presumably give nodal errors of that magnitude. Therefore while QMC 
calculation of even larger systems are possible as the systems size increases so does the 
nodal error. Methods to improve the nodal structure are required to achieve a total energy 



12 



resolution of 0.1 eV. We believe the above observations will guide method selection for exci- 
tations in nanostructures. Clearly, the robust design of optical nanoarchitectures will require 
further developments in many-body theories such as GW-BSE and QMC. We also encour- 
age more precise experimental determination of electronic excitations in well-characterized 
nanostructures to provide a more severe test to theory. 
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FIG. 1: (Color online) Structures of the fullerenes with corresponding point group representations 
in parentheses: C 2 o (I/i), C 2 4 (D 6 d), C 50 (D 5ft ), C 60 (I/O, C 70 (D 5h ), C 8 o (I/i), C 80 (D 5d ). 
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